Instability of vortex array and transitions to turbulent states in rotating helium II 
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We consider superfluid helium inside a container which rotates at constant angular velocity and 
investigate numerically the stability of the array of quantized vortices in the presence of an imposed 
axial counterflow. This problem was studied experimentally by Swanson et at, who reported evi- 
dence of instabilities at increasing axial flow but were not able to explain their nature. We find that 
Kelvin waves on individual vortices become unstable and grow in amplitude, until the amplitude 
of the waves becomes large enough that vortex reconnections take place and the vortex array is 
destabilized. The eventual nonlinear saturation of the instability consists of a turbulent tangle of 
quantized vortices which is strongly polarized. The computed results compare well with the exper- 
iments. Finally we suggest a theoretical explanation for the second instability which was observed 
at higher values of the axial flow. 

PACS numbers: PACS 67.40.Vs, 47.37.+q, 03.75.Lm 



I. INTRODUCTION 

The work described in this article is concerned with the 
stability of a superfluid vortex array. It is well known 12 
that, if helium II is rotated at constant angular velocity 
f2, an array of superfluid vortex lines is created. The 
vortices are aligned along the axis of rotation and form 
an array with areal density given by 
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where n = h/m = 9.97 x 10~ 4 cm 2 /sec is the quantum 
of circulation, h is Plank's constant and m the mass of 
one helium atom. Equation (JTJ is valid provided that f2 
exceeds a small critical valued. Rotation frequencies of 
the order of 1 Hz are easily achieved in a laboratory, and 
correspond to real density of the order of 10 3 cm~ 2 . 

It is also well known that a superfluid vortex line be- 
comes unstable in the presence of normal fluid in the di- 
rection parallel to the axis of the vortex. This instability, 
hereafter referred to as the Donnelly- Glaberson (DG) in- 
stability, was first observed experimentally by Cheng et 
al~ and then explained by Glaberson et al^. Physically, 
the DG instability takes the form of Kelvin waves (helical 
displacements of the vortex core) which grow exponen- 
tially with time. 

In this paper we use an imposed axial flow to trigger 
the DG instability and study the transition from order to 
disorder in an array of quantized vortex lines. It is useful 
to remark here that, since the growth of Kelvin waves 
takes place at the expense of normal fluid's energy, un- 
derstanding the DG instability is also relevant to the 
balance of energy between normal fluid and superfluid in 
helium II turbulence, a problem which is attracting cur- 
rent experimenta£SSiS and theoreticaliiASA 3 ^ attention. 



The article is organized in the following way. In sec- 
tion 2 we describe the rotating counterflow configura- 
tion, which is relevant to both theory and experiment. 
In section 3 we summarize experimental results obtained 
by Swanson, Barenghi and Donnelly 15 . They discovered 
that the DG mechanism can destabilize the superfluid 
vortex array and revealed the existence of two different 
superfluid states at increasing values of the driving ax- 
ial flow beyond the DG transition. Until now, the ac- 
tual physical nature of these two states has been a mys- 
tery, and it is the aim of our work to shed light into 
this problem. In section 4 we set up the formulation of 
vortex dynamics in the rotating frame which generalizes 
the previous approach of Schwarai& and which we use in 
our numerical calculations. Section 5 is devoted to the 
DG instability. What happens beyond the DG instability 
cannot be predicted by linear stability theory and must 
be investigated by direct nonlinear computation, which 
is done in section 6. In section 7 we tackle the transition 
to the second turbulent state observed by Swanson et al. 
Finally, section 8 draws the conclusions. 



II. ROTATING COUNTERFLOW 

In order to study the stability of the rotating super- 
fluid vortex array, we consider the configuration which 
is schematically shown in Fig. 1. A channel, which is 
closed at one end and open to the helium bath at the 
other end, is placed on a table which can be rotated at 
an assigned angular velocity fl. At the closed end of the 
channel a resistor dissipates a known heat flux Q. 

First let us consider what happens in the absence of 
rotation (Q — 0). Since only the normal fluid carries 
entropy, then Q — pTSV n , where T is the absolute tem- 
perature, S the specific entropy, p = p s + p n the total 
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at given values of f2 and V ns , and, if the array is unsta- 
ble, of whether the vortex line density L is the sum of 
Eq. and Eq. or not. Their experimental results 
are described in the next section. 

It is important to remark that, in principle, one can 
also study the stability of a vortex array in the presence of 
a mass flow rather than of a heat current. Similarly, one 
can study the effects of rotation upon the turbulence of 
helium II created by towing a grid or rotating a propellers 
rather than upon counterflow turbulence. The reason for 
which we have chosen to restrict our investigation to the 
case of a heat current is twofold: firstly, the experimental 
data of Swanson et ate are available; secondly, at least at 
small heat currents^, the turbulent superfluid tangle is 
homogeneous and almost isotropic and we do not have to 
worry about large scale motion and eddies of the normal 
fluid. 



FIG. 1: Schematic rotating counterflow apparatus. 

density of helium II, p s the superfluid density and p n 
the normal fluid density. We call V n and V s respectively 
the normal fluid and the superfluid velocity fields in the 
direction along the channel, averaged over the channel's 
cross section. The total mass flux p s V s + p n V n is zero 
because one end of the channel is closed. The resulting 
counterflow velocity V ns = V n — V s which is induced along 
the channel is therefore proportional to the applied heat 
flux: 



It is known from experimentsi£A& and numerical 
simulations^ that, if Q (hence V ns ) exceeds a critical 
value, a turbulent tangle of quantized vortex lines is cre- 
ated. The tangle is homogeneous and isotropic (neglect- 
ing a small degree of anisotropy induced by the direction 
of the imposed hear current). The intensity of the tur- 
bulence is measured by the vortex line density (length of 
vortex line per unit volume) which is experimentally de- 
termined by monitoring the extra attenuation of second 
sound. It is found that the vortex line density has the 
form 

iflow = JhV^ s , (3) 

where •jh is a temperature dependent coefficienlte. 

Let us consider now the case in which the heat flux is 
applied in the rotating frame (fi ^ 0). We have now 
two effects which compete with each other: rotation, 
which favours the creation of an ordered array of vor- 
tices aligned along the direction of the axis of rotation, 
and counterflow, which favours the creation of a disor- 
dered tangle. Swanson et aliU were the first to address 
the problem of whether the vortex array is stable or not 



III. THE EXPERIMENT 

The rotating counterflow apparatus of Swanson et ate 
consisted of a 40 cm long vertical channel with lxl cm 2 
square cross section. At the closed end (as shown in Fig. 
1) a resistor dissipated a known heat flux Q and induced 
relative motion V ns of the two fluid components. The 
vortex line density L was measured by pairs of second 
sound transducers located along the channel. The entire 
apparatus was set up on a rotating cryostat, so that it 
was possible to create vortex lines by either rotation or 
counterflow, or by any combination of them. The vortex 
line density was calculated from a measurement of the at- 
tenuation of second sound resonances and its calibration 
against the known density in rotation^. 

The experiment was performed at T — 1.65 K. In the 
presence of both rotation and counterflow three distinct 
flow states were observed, as shown in Fig. 2. The three 
states are separated by two critical counterflow veloci- 
ties V c i and V C 2, which are respectively the boundaries 
between the primary state and the secondary state, and 
between the secondary state and the tertiary state. The 
results of the experiment can be summarized as it follows: 

• Primary state 

In the first region of Figure 2 at the left of V c \ the vor- 
tex line density is independent of the small values of V ns 
involved and agrees with Eq. ■ This region clearly cor- 
responds to an ordered vortex array, and the counterflow 
current V ns is not strong enough to destabilize it. 

• Transition from primary state to secondary 
state 

Swanson et aL 15 noticed that the values of the first 
critical velocity V c i are consistent with the DG instabil- 
ity. This means that at V ns — Vbg the axial flow is so 
strong that Kelvin waves of infinitesimal amplitude be- 
come unstable. 

• Secondary state 

Because of the lack of direct flow visualization in he- 
lium II, the nature of the flow past the instability (V ns > 
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Vci) was not clear to Swanson et a&^. The only infor- 
mation which they could recover by the second sound 
measurement technique was that rotation added fewer 
than the expected 2Q / k vortex lines to those made by 
the counterflow current. 

• Transition from secondary state to tertiary state 

The existence of a second critical velocity V C 2 was un- 
expected. The nature of the transition at V ns = V C 2 and 
which kind of flow exists in the third region (V ns > V C 2) 
were a mystery. 
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FIG. 2: Vortex line density L observed by Swanson et alw^ 
in the presence of a counterflow heat current V„ s at various 
rotation rates at T = 1.65. The solid lines represent fits to the 
two observed critical velocities V c i and V C 2- The experimental 
uncertainties are about 1/3 of the symbol size. 



IV. VORTEX DYNAMICS IN A ROTATING 
FRAME 

The vortex filament model is very useful to study the 
motion of superfluid 4 He because the vortex core ra- 
dius do ~ 10~ 8 cm is microscopic, hence much smaller 
than any flow scales of interest. Moreover, unlike what 
happens in classical fluid dynamics, the circulation k = 
9.97x 10 -4 cm 2 /sec is fixed by quantum constraint, which 
simplifies the model ever further. 

Helmholtz's theorem for a perfect fluid states that a 
vortex moves with self - induced velocity at each place 
produced by the shape of the vortex itself. Therefore 
the velocity s*o of a vortex filament at the point s in the 
absence of mutual friction is governed by the Biot-Savart 
law and can be expressed asi& 
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Here the vortex filament is represented by the parametric 
equation s — s(£,i). The first term means the localized 
induction velocity, where the symbols 1+ and £_ are the 
lengths of the two line elements which are adjacent to a 
given point after discretization of the filament, and the 
prime denotes differentiation of s with respect to the arc 
length £. The second term represents the nonlocal field 
obtained by carrying out the integral along the rest of 
the filament on which s± refers to a point. 

If the temperature is finite, the normal fluid fraction is 
non-zero and its effects must be taken into account. The 
normal fluid induces a mutual friction force which drags 
the vortex core of a superfluid vortex filament for which 
the velocity of point s is given by 



s + as x (v ns - s ) - as x [s x (v ns - s )], (5) 



where a and a' are known temperature-dependent fric- 
tion coefficients and so is calculated from Eq. More 
details of the numerical method and how it is imple- 
mented are described in RefiSi. 

In order to make progress in our problem, we need to 
generalize this vortex dynamics approach to a rotating 
frame 2 ^. The natural way to perform the calculation in 
a rotating frame would require to consider a cylindrical 
container. We do not follow this approach for two rea- 
sons. Firstly, our formulation is implemented using the 
full Biot - Savart law, not the localized-induction approx- 
imation often used in the literature. This would require 
to place image vortices beyond the solid boundary to im- 
pose the condition of no flow across it. This is easily done 
in cartesian (cubic) geometry, but it is difficult to do in 
cylindrical geometry, Secondly, the original experiment 
by Swanson et ali^ was carried out in a rotating channel 
with a square cross section. 

In a rotating vessel the equation of motion of vortices 
is modified by two effects. The first is the force acting 
upon the vortex due to the rotation. According to the 
Helmholtz's theorem, the generalized force acting upon 
the vortex is balanced by the Magnus force: 
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where F' = F — £1 • M is the free energy of a system 
in a frame rotating around a fixed axis with the angular 
velocity Q and the angular momentum M. Taking the 
vector product of Eq.JfJJ with s 1 , we obtain the velocity 
s*o • The first term F due to the kinetic energy of vortices 
gives that Biot - Savart law, and the second term Q, ■ M 
leads to the velocity s ro t of the vortex caused by the 
rotation: 
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with B — r—s. The second effect is the supcrffow induced 
by the rotating vessel. For a perfect fluid we know the 
analytical solution of the velocity inside a cube of size D 
rotating with the angular velocity = QaQ. 
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with X = (2n + l)irx/D and V = (2n + l)nY/D. In a 
rotating frame these terms are added to the velocity s*o 
without the mutual friction, so Eq. {3} is replaced by 

+ Srot+Ucub- (10) 

Some important quantities useful for characterizing the 
rotating tangle will be introduced. The vortex line den- 
sity is 



(11) 



where the integral is made along all vortices in the sam- 
ple volume A. The polarization of the tangle may be 
measured by the quantity 



<4>=-L / ,/^(<|.-.. 



(12) 



as a function of time. 

The actual numerical technique used to perform the 
simulation has been alreay described—. Here it is enough 
to say that a vortex filament is represented by a sin- 
gle string of points at a distance A£ apart. When 
two vortices approach within A£, they are assumed to 
reconnects^. The computational sample is taken to be 
a cube of size D — 1.0 cm. We adopt periodic bound- 
ary conditions along the rotating axis and rigid bound- 
ary conditions at the side walls. All calculations are 



made under the fully Biot-Savart law, placing image vor- 
tices beyond the solid boundaries. The space resolu- 
tion is A£ = 1.83 x 10~ 2 cm and the time resolution 
is At = 4.0 x 10~ 3 sec. All results presented in this pa- 
per refer to calculations made in the rotating frame. To 
make comparison with the experiment^, we use a = 0.1 
and a' = at the temperature T = 1.6K. The uniform 
counterflow v ns is applied along the z axis. 



V. THE DG INSTABILITY 

Swanson et al. 15 found that the first critical velocity 
V c i was proportional to fi 1 / 2 ; this functional dependence 
and the actual numerical values were consistent with in- 
terpreting the transition at V ns — V c \ as the DG insta- 
bility of Kelvin waves. Glaberson et al~ considered an 
array of quantized vortices (which they modelled as a 
continuum) inside a container rotating at an angular ve- 
locity ft. They found that, in the absence of friction, the 
dispersion relation of a Kelvin wave of wavenumber k is 



uj = 2ft + uk 2 



(13) 



where uj is the angular frequency of the Kelvin wave, v 
is given by 



fln(^), 

47T do 



(14) 



where b ~ L^ 1 / 2 is the average distance between vortices. 

Glaberson et aL 5 showed that the dispersion law l|13|) 
has a critical velocity 



Vdg = 7 

K rain 

at the critical wavenumber 
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If the axial flow V ns exceeds Vdg for some value of k, 
then Kelvin waves with that wavenumber k (which are 
always present at very small amplitude due to thermal 
excitations and mechanical vibrations) become unstable 
and grow exponentially in time. Physically, the phase 
velocity of the mode k is equal to the axial flow, so energy 
is fed into the Kelvin wave by the normal flow. 

Figure 3 illustrates the DG instability. The compu- 
tations were performed in a periodic box of size 1 cm 
in a reference frame rotating with angular velocity Q — 
9.97 x 10~ 3 rad/sec, for which V DG = 0.010 cm/sec. Fig- 
ure 3a confirms that when V ns — 0.008 cm/sec < Vdg 
the vortex lines remain stable. Figure 3b shows that, at 
V ns = 0.015 cm/sec > Vdg, Kelvin waves become un- 
stable and grow, as predicted. Figures 3c,d,e and f show 
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FIG. 3: Numerical simulations of the Donnelly-Glaberson in- 
stability at tt = 9.97 x 10~ 3 rad/sec, T = 1.6K. Snapshots 
of vortex configurations at the following counterflow veloci- 
ties V ns ' (a): Vns = 0.008cm/sec; (b): V ns = 0.015cm/sec; 
(c): Vns = 0.03cm/sec; (d): V„ s = 0.05cm/sec; (e): Vns = 
0.06cm/sec; (f): V ns = 0.08cm/sec. 

that Kelvin waves of larger wavenumber become unstable 
at higher counterflow velocity. 

Linear stability theory^ can only predict two quanti- 
ties: the first is the critical value of the driving parameter 
{Vdg m our case) at which a given state (the vortex ar- 
ray in our case) becomes unstable because infinitesimal 
perturbations grow rather than decay; the second is the 
exponential growth or decay rate of these perturbations 
for a given value of the driving parameter. Therefore the 
linear stability theory of Glaberson^ cannot answer the 
question of what is the new solution which grows beyond 
the DG instability: to determine this new solution (what 
we call the secondary state in section 3) we must solve the 
governing nonlinear equations of motion, which is what 
we do in the next section. 



VI. ROTATING TURBULENCE 

Because of the computational cost of the Biot-Savart 
law, it is not practically possible to compute vortex tan- 
gle with densities which are as high (L = 0(lO 4 )cm~ 2 ) as 



those achieved in the experiment. Nevertheless, numer- 
ical simulations performed at smaller, computationally 
realistic values of L are sufficient to shed light into the 
physical processes involved. Some results which we de- 
scribe have been already presented in preliminary formSi; 
together with more recent computer simulations, the pic- 
ture which emerges and which we present here gives a 
good understanding of the experimental findings of Swan- 
son et ate, at least as far as the transition to the sec- 
ondary flow and the secondary flow itself are concerned. 

The time sequence contained in Fig. 4 illustrates the 
evolution of a vortex array at a relatively small angular 
velocity f2 = 9.97 x 10 _3 rad/sec, in the presence of the 
counterflow V ns = 0.08cm/sec. Figure 4a shows the ini- 
tial N = 8 parallel vortex lines at t = 0. The vortices 
have been seeded with small random perturbations to 
make the simulation more realistic. The absence of these 
perturbations would make the phase of the Kelvin waves 
synchronize on all vortices to delay reconnections. As the 
evolution proceeds, perturbations with high wavenum- 
bers are damped by the mutual friction, whereas per- 
turbations which are linearly DG-unstable grow expo- 
nentially, hence Kelvin waves become visible (Fig. 4b). 
When the amplitude of the Kelvin waves becomes of the 
order of the average vortex separation, reconnections take 
place (Fig. 4c). The resulting vortex loops disturb the 
initial vortex array, leading to an apparently random vor- 
tex tangle (Fig. 4d). After the initial exponential growth 
(which is predicted by the theory of the DG instabil- 
ity), nonlinear effects (vortex interactions and vortex re- 
connections) become important and nonlinear saturation 
takes place. 

Figure 5 shows a similar time sequence at the same 
counterflow velocity V ns = 0.08cm/sec but at higher ro- 
tation rate = 4.98 x 10 _2 rad/sec. In this case we have 
N = 33 initial parallel vortices (Fig. 5a). At t — 12sec 
(Fig. 5b) it is still N — 33. Then the amplitude of the 
Kelvin waves becomes so large that lots of reconnections 
take place and N increases; for example, we have N = 83 
at t = 160sec (Fig. 5f). 

It is instructive to compare these results with ordinary 
counterflow in the absence of rotation. Figure 6 shows a 
vortex tangle obtained for = and V ns — 0.08cm/sec. 
The dynamics starts from N = 6 vortex rings. It has 
been known since the early work of SchwaraAS that the 
resulting tangle does not depend on the initial condition. 
In this particular simulation the vortices develop to a 
turbulent tangle. 

Figure 7 shows that in all three cases (small rotation, 
large rotation, no rotation) the vortices, after an ini- 
tial transient, saturate to a statistically steady, turbulent 
state, which is characterized by a certain average value 
of L. In the case of f2 ^ (Fig. 6a and b), it is apparent 
that the initial growth is exponential, which confirms the 
occurrence of a linear instability. 

Looking carefully at the saturated tangle at higher ro- 
tation (Fig. 5f) we notice that there are more loops ori- 
ented vertically than horizontally. The effect is not visi- 



(a) (b) 



(c) (d) 




FIG. 4: Numerical simulation of rotating turbulence at T = 
1.6K, n = 9.97 x 10" 3 rad/sec and V ns = 0.08cm/sec. Com- 
puted vortex tangle at the following times: (a): t=0rm sec; 
(b): t=16sec; (c): t=28sec; (d): t=36sec; (e): t=80sec; (f): 
t=600sec. 



(a) (b) 




(c) (d) 




(e) (f) 

FIG. 5: Numerical simulation of rotating turbulence at T = 
1.6K, n = 4.98 x 10" 2 rad/sec and V ns = 0.08cm/sec. Com- 
puted vortex tangle at the following times: (a): t=0rm sec; 
(b): t=12sec; (c): t=20sec; (d): t=28sec; (e): t=40sec; (f): 
t=160sec. 



ble at lower rotation (Fig. 4f) and at zero rotation (Fig. 
6f). The degree of polarization of the tangle is repre- 
sented by < s' z > of Eq. i|12fl . This quantity captures the 
difference between a vortex array (for which < s' z >= 1 
because all lines are aligned in the +z direction) and 
a random vortex tangle (for which < s' z >= because 
there is an equal amount of vorticity in the +z and — z 
directions). Figure 8 shows how < s' z > changes with 
time in the three cases (small rotation, large rotation, no 
rotation) considered. The quantities of interest are the 
values of < s' z > at large times in the saturated regimes. 
In the absence of rotation (Fig. 8c) < s' z > is negligible 
but not exactly zero (< s' z >~ 0.01), certainly because 
the driving counterflow is along the z direction. This 
small anisotropy of the counterflow tangle has been al- 
ready reported in the literature^. At small rotation (Fig. 
8a) there is a small but finite polarization < s' z >« 0.15, 
whereas at higher rotation (Fig. 8b) the polarization is 
significant (< s' z >m 0.45) - in fact it is even visible with 
the naked eye (Fig. 5f). 

Figure 9 shows the calculated dependence of the vortex 
line density L on the counterflow velocity V ns at differ- 
ent rotation rates fl. The figure shows a dependence of 
L on V ns which is similar to what appears in the Fig. 1 
of the paper by Swanson et a! 15 . The only difference is 
that the scale of the axes in the paper by Swanson al. is 



bigger - in this particular figure they report vortex line 
densities as high as L w 2500cm -2 , whereas our calcu- 
lations are limited to L w 80cm -2 . Despite the lack of 
overlap between the experimental and numerical ranges, 
there is clear qualitative similarity between the figures. 
It is apparent that the critical velocity beyond which L 
increases with V ns is much reduced by the presence of 
rotation, which is consistent with the observations. 

Figure 10 shows the calculated polarization < s' z > 
as a function of counterflow velocity V ns at different ro- 
tation rates O. It is apparent that the polarization de- 
creases with the counterflow velocity and increases with 
the rotation, which shows the competition between order 
induced by rotation and disorder induced by flow. 

We conclude that the nonlinear saturation which takes 
place beyond the DG instability and which was observed 
by Swanson et al*^ is a state of "polarized" turbulence. 



VII. THE SECOND CRITICAL VELOCITY 

In this section we propose a qualitative theory for the 
second critical velocity which was observed, but not ex- 
plained, by Swanson et al. 15 . Unfortunately, this region 
of parameter space cannot be explored directly using nu- 
merical methods, due to the larger vortex line densities. 
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(e) 



(f) 



FIG. 6: Numerical simulation of counterflow turbulence at 
T — 1.6K in the absence of rotation (fi = 0) for V na = 
0.08cm/scc. Computed vortex tangle at the following times: 
(a): t=0rm sec; (b): t=120sec; (c): t=360sec; (d): t=520sec; 
(e): t=680sec; (f): t=1160sec. 




FIG. 7: Vortex line density L vs time t at T = 1.6K and 
V ns = 0.08cm/sec for: (a): Q = 9.97 x 10 _3 rad/sec; (b): 
fi = 4.98 x !CT 2 rad/sec; (c): Q = Orad/sec. 




FIG. 8: Tangle's polarization < s' z > vs time t at T = 1.6K 
and V ns = 0.08cm/sec for: (a): fl = 9.97 x l(T 3 rad/scc; (b): 
n = 4.98 x l(T 2 rad/sec. (c): Q = Orad/sec. 
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FIG. 9: Vortex line density L vs V^ s at T = 1.6K for 
fi = (white square), 57 = 9.97 x 10~ 3 rad/sec (black square), 
n = 2.99 x l(T 2 rad/sec (circle) and Q = 4.98 x 10 _2 rad/sec 
(triangle) . The dotted lines are guides to the eye. 



We consider the following idealized model of the region 
V ns > V C 2 ■ We represent the polarized tangle as the com- 
bination of a vortex array and a number of vortex loops 
which are randomly oriented, so that the combined sys- 
tem is in balance and has the necessary amount of length 
and polarization. Let T\ be the characteristic timescale of 
the growing Kelvin waves which are induced on the vor- 
tex array by the DG instability. The typical lifetime T2 of 
the vortex loops will be determined by the friction with 
the normal fluid and by the relative orientation with re- 
spect to the counterflow. If r 2 > T\ then the vortex loops 
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FIG. 10: Tangle's polarization < s' z > vs V% a at T — 1.6K 
for n = 9.97 x 10~ 3 rad/sec (square), fi = 2.99 x 10~ 2 rad/sec 
(circle) and = 4.98 x 10 _2 rad/sec (triangle). The dotted 
lines are guides to the eye. 



will not have enough time to shrink significantly before 
more loops are introduced by vortex reconnections in- 
duced by growing Kelvin waves. The combined system 
consisting of the vortex array and the vortex loops will 
not be balanced any longer because randomness is intro- 
duced by vortex reconnections at a rate which is faster 
than the rate at which loops are destroyed by friction. In 
conclusion, we expect that the vortex tangle will be tur- 
bulent and unpolarized. According to this scenario, the 
order of magnitude of the critical velocity V c i is given by 
the condition 



n = r 2 . 



(17) 



First we estimate t\ using a simple model. For the 
sake of simplicity we assume an isolated vortex line of 
helical shape s = (e cos 4>; e sin <j>; z) where cj> — kz—tut and 
e « 1, hence z ~ £ is the arc- length. The tangent unit 
vector is s* = ds/d£ ~ ds/dz = (— fee sin 0; kecostp; 1) 
and s" = (— fc 2 ecos</>; — fc 2 esin</>;0). Using the localized- 
induction approximation, the self-induced velocity of the 
line at the point s is given by 



«ns = (0;0;V^ S ) and e = e(i), we find that de/dt — 
a(kV ns — v'k 2 )e hence e(i) = e(0)exp(at) where the 
growth rate is a = a(kV ns — v'k 2 ). Given V ns , the 
largest growth rate occurs for k = V ns /(2v') and takes 
a = orV^/(4i/), for which we conclude that 



1 W 

n = — 

a 



(19) 



To estimate t 2 we approximate the vortex loops as vor- 
tex rings of radius approximately determined by the aver- 
age vortex spacing S w £ -1 ' 2 . The characteristic lifetime 
of a ring of radius R in the presence of friction is 2 - 



Tl = 



2p s irR 2 

7^2 



(20) 



where £2 = In [(8-R/ao) — 1/2] and 7 is a known friction 
cocfficien1j2£. Setting 2R = S = L -1 / 2 , we conclude that 
the polarized tangle is unstable if 



L < C 2 VJ 



with 



C 2 



2 7 r£ 1 £ 2 ' 



(21) 



(22) 



Equation l|21() has the same dependence of L on V ns as 
that observed experimentally. At T — 1.65K we have2£ 
p s = 0.1168g/cm 3 , 7 = 1.3 x l(T 5 g/cm sec, a = 0.11. 
Since clq k 10 _8 cm and the slowly varying logarithm 
terms are approximately C\ « £2 ~ 10, we conclude 
that C2 ~ 5 x 10 4 cm~ 4 sec 2 , which is of the same order 
of magnitude of the value C 2 = 16 x 10 4 cm" 
by Swanson et ai&. 



Vi = v s x s , (18) 

where v 1 = kCi/(4lTt) and the slowly varying term C\ 
!*i [l/(fcao)] is assumed constant. Neglecting higher order 
terms in e we have Vi = j/'fc 2 e(sin 0; — cos </>; 0). 

In the absence of friction the equation of motion is 
simply ds/dt = hence, assuming that e is constant, 
we find that the Kelvin wave oscillates with angular fre- 
quency uj — v'k 2 . This result differs from Glaberson's 
Eq. i|13|) because we perturbed a single vortex line rather 
than a continuum of vorticity 2f2 described by the Hall- 
Vinen equations in the rotating frame (hence the pres- 
ence of a different upper cutoff which makes v' different 
from v and the contribution 2f2 to to). 

In the presence of friction, neglecting the small mu- 
tual friction coefficient a 1 for simplicity, the equation 
of motion is ds/dt = + as 1 x (v ns — v^). Assuming 



VIII. CONCLUSIONS 

In conclusion, we have studied the stability of a super- 
fluid vortex array in the presence of an applied coun- 
tcrflow, giving answers to some questions which were 
first asked by the pioneering experiment of Swanson et 
(ili^L. After investigating the DG instability, we have de- 
termined the existence of a new state of superfluid tur- 
bulence (polarized turbulence) which is characterized by 
two statistically steady state properties, the vortex line 
density and the degree of polarization. Although our 
computed range of vortex line densities does not overlap 
with the much higher values obtained in the experiment, 
we find the same qualitative dependence of vortex line 
density versus counterflow velocity at different rotations. 
We have also made some qualitative progress to under- 
stand what happens beyond U c2 . 
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Further work with more computing power will hope- 
fully investigate other aspects of the problem, particu- 
larly what happens at high counterflow velocities and 
line densities. We also hope that our work will stimulate 
more experiments on this problem. 

It is somewhat surprising that so little is known about 
the destabilization of a rotational vortex array by an im- 
posed counterflow. For example, it should be possible to 
observe the polarization of turbulence by using simulta- 
neous measurements of second sound attenuation along 
and across the rotation axis. 

Finally, our work should be of interest to other investi- 
gations of vortex arrays and how they can be destabilized 
in other systems, ranging from superfluid 3 He2£to atomic 
Bose-Einstein condensates- 30 . It is also worth noticing 
that this study has revealed the crossover of the dimen- 



sionality of vortex systems. If one considers the three 
regimes in Fig. 2 one notices that, at a fixed value of 
V ns , increasing the rotation rate makes the vortices po- 
larized, changing the dynamics from three-dimensional 
to two-dimensional. This reduction of the dimension- 
ality of turbulence has been observed in classical fluid 
mechanics^!. 
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